Gravity in the brane-world for two-branes model with stabihzed modulus 
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We present a complete scheme to discuss linear perturbations in the two-branes model of the 
Randall and Sundrum scenario with the stabilization mechanism proposed by Goldberger and Wise. 
We confirm that under the approximation of zero mode truncation the induced metric on the branes 
reproduces that of the usual 4-dimensional Einstein gravity. We also present formulas to evaluate the 
mass spectrum and the contribution to the metric perturbations from all the Kaluza-Klein modes. 
We also conjecture that the model has tachyonic modes unless the background configuration for the 
bulk scalar field introduced to stabilize the distance between the two branes is monotonic in the 
fifth dimension. 
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I. INTRODUCTION 

Recently, there has been a growing interest in consider- 
ing extra-dimensions in non-trivial form One im- 
portant suggestion was done by Randall and Sundrum 
1^. They realized that if one considers 5-dimensional 
gravity with a negative cosmological constant bounded 
by two positive and negative tension branes, the induced 
gravity on the negative tension brane can possess a very 
small gravitational constant as compared with the en- 
ergy scale introduced in the original Lagrangian. This 
fact potentially gives a solution to the hierarchy problem 
of explaining the extraordinary weakness of the gravita- 
tional coupling. 

Later, the same authors pointed out that the model 
with the positive tension brane alone is also possible . 
In this case, the extension of the extra dimension is infi- 
nite. Nevertheless, the induced gravity on the brane can 
be expected to mimic Einstein gravity. 

There are many discussions about the cosmology based 
on these scenarios The behavior of gravity on 

these models has been also investigated by many authors 
pl|-|r7t . In this direction, an explicit method to deal with 
the linear order metric perturbations induced by the mat- 
ter fields confined on the branes was developed in the pa- 
per by Garriga and Tanaka (Paper I). In that paper, 
considering the zero mode truncation approximation, it 
was shown that the induced gravity on the branes be- 
comes of the Brans-Dicke type in the two-branes model. 
On the other hand, Einstein gravity is recovered in the 
single-brane model. 

However, the analysis of the two-branes model in Pa- 
per I was not complete. The stabilization mechanism of 
the distance between the two branes by means of 

a bulk scalar field was not taken into account. Moreover, 
the Kaluza-Klein contribution was not estimated for the 
two-branes model. In this paper we consider a model 
with stabilization mechanism, and we confirm that the 4- 
dimensional Einstein gravity is recovered on both branes 



under the approximation of zero mode truncation. An 
approximate formula for the mass of the lowest massive 
mode in the scalar-type perturbations is also obtained. 
Furthermore, we present a method to evaluate the met- 
ric perturbations taking into account all the KK contri- 
butions under a contact interaction approximation. This 
approximation is valid if the source which excites the KK 
modes is smoothly distributed compared with the scale 
corresponding to the mass of the lowest massive mode. 

This paper is organized as follows. In Sec. II we de- 
scribe the background model that we consider, and derive 
the basic equations for the linear perturbations on it. In 
Sec. Ill we discuss the mass spectrum of the perturba- 
tions. We identify the mode functions for the massless 
degrees of freedom, and also derive the formula for the 
mass and the mode function of the lowest massive mode. 
In Sec. IV the mechanism to recover Einstein gravity is 
explained. This result is expected from the notion of the 
mass spectrum but its derivation is not so trivial when 
we consider the explicit construction of metric perturba- 
tions. We shall find that the contribution from the mas- 
sive modes of the scalar-type perturbations must be also 
taken into account in some sense. For this purpose, we 
introduce a contact interaction approximation. In Sec. V 
we apply the same technique to the KK modes of the 
tensor-type perturbations to complete the description of 
the linear perturbations of this system. Section VI is 
devoted to summary. 



II. MODEL AND BASIC EQUATIONS 

We consider linear perturbations of the model pro- 
posed by Randall and Sundrum |^ with the stabiliza- 
tion mechanism by Goldberger and Wise ||l8| . The fields 
existing on the 5-dimensional spacetime (bulk) are the 
usual 5-dimensional gravity with a negative cosmological 
constant A, and a 5-dimensional scalar field f introduced 
to stabilize the distance between the two branes. We de- 



note the 5-dimensional gravitational constant by G^. The 
unperturbed metric is supposed to take the form, 



a^iy)r]fj,^dx'^dx'' + dy^ 



(2.1) 



where rjf^i, is the 4-dimensional Minkowski metric with 
( — h ++) signature. We use = a?{y)'q^i, to raise 4- 
dimensional tensor indices. The y-direction is bounded 
by two branes located at y — y*-^-*. On these two branes, 
Z2-symmetry is imposed. The Lagrangian for the bulk 
scalar field is 



V-bM- E^*'^M%-y^"^)• 
(2.2) 



<T = ± 



For most of the present analysis, we do not need to spec- 
ify the explicit form of the potentials Vsif ) and F^^^ {(p). 
In the bulk, the background scale factor and the scalar 
field, {a,(po), must satisfy 

Hiy) = -^^liy), 



HHy) 



6 V2' 



ipl{y)~VBiMv))- ^^ 'A 



(2.3) 



where H{y) := d{y)/a{y) « — \/— A/6 and k = SttGs. 
Ordinary matter fields reside on both branes, and the val- 
ues of the 4-dimensional vacuum energy on both branes 
are adjusted to realize a static background configuration. 
This way of model construction is different from the stan- 
dard approach jl^ which assumes the model potential 
and the values of the vacuum energy on the branes from 
the very beginning. In the present approach, we obtain 
a constraint on the model potential by fixing the dis- 
tance between the two branes. Conversely, in the stan- 
dard approach, the stabilized distance is determined for 
each given model (although some tuning of one of the 
model parameters is necessary to realize a static config- 
uration). We call the brane &i y — y^^\y''^^) the posi- 
tive (negative) tension brane, and adopt the convention 

y( + )<y(-). 

We first consider the metric perturbation 5{ds'^) — 
habdx"" dx^ and the scalar field perturbation S^p in the 
bulk. We are using the convention that Latin indices 
run through 0, 1, 2, 3, 5, while Greek indices run through 
0, 1, 2, 3. In the bulk, we can always impose the "Newton 
gauge" condition. 



h55 = 20™ , /l5^ = 0, 



5ip 



[dy + 2H] 



(2.4) 



where h'l^^^ satisfies the transverse-traceless condition, 
and (/)^ is the 5-dimensional "Newtonian potential" . 
The equation for /i^^ in the bulk is given by 



1 



□ (4) ^ liTT) 



hZ^^ = 0, 



(2.5) 



where is the 4-dimensional d'Alembertian operator 
with respect to ry^^ , and we have introduced the operator 



For the scalar-type perturbations, we obtain 



l(4) 



L 



= 0, 



(2.6) 



(2.7) 



where L^"^ ^ is the operator defined by 



"^Q^v-TT^^ya - —a ipo- 

(Jj Y^n O 



(2., 



As discussed in Paper I, the y = constant hypersurface 
does not correspond to the location of the branes in this 
gauge. Following it, we introduce two Gaussian normal 
coordinate systems near both branes, respectively. We 
denote quantities in these coordinate systems by associ- 
ating a bar, such as y. 

Then, the junction condition for the bulk scalar field 
is given by ±2(^ = y(^^((^) at y = y'-^-'i. Hence, its 
perturbation become 

±25(p^V"^^\ipa)5(p, (y = y(±)±). (2.9) 

The junction condition for the metric perturbation is 



2k 



(±) 



T^7^.</'o5^, (y = y(±)±). (2.10) 



Here, we have introduced the energy momentum tensor 
of the matter fields confined on the branes, Tjhj\ and 
rp{±) jf^'^T^p, satisfies the 4-dimensional con- 

servation law r^jy''' = 0. Notice that there appears a 
contribution from the perturbation of the scalar field f 
which was not present in the models discussed in Paper I. 

Whereas the junction conditions have been easily de- 
rived using Gaussian normal coordinates, the equations 
of motion for the perturbations are simpler in the "New- 
ton gauge" . So we need to consider the gauge transforma- 
tion which relates normal coordinates and the "Newton 
gauge" . This is given by hab = hab + S,a-b + S.b;a with 



%^{y')dy'^ r ^""{yW + iti). 
\^''{y')dy' r <l>l{y")dy" 



(±) 



±)' 



(2.11) 
(2.12) 



2 



where ^^^^^ and are independent of y. Then, we have 



d(p{y) = Sip{y) - ipoiy) 



(±) 



dy' 



4>Uy")dy" 



(2.13) 



Here we should stress that the arguments in the Lh.s. are 
not y but y. 



Combining Eqs.( 2.10 ) and ( 2.13 ), the junction condi- 
tion in the "Newton gauge" for the TT part is obtained 
as 



±(a,-2i/)e/^ = -KSW, (y = 2;W±) 



where we have defined 



(2.14) 



T, 



1 



(±) 



1 



(±),P 



(2.15) 



tion for }ij[J"\ 



Combining Eq. (2.14) with Eq. ( p.5[ ), we obtain the equa- 
on fc 

□ (4) 



/zf/) = -2.^I](^)<5(y-j/('^)). (2.16) 



(T = ± 



(TT) 

The resulting /ij^^ is automatically consistent with 
the traceless condition, but the transverse condition gives 
us the equation which determines 



(±)' 



-□(4)^5 



(±) 



(2.17) 



wher e we have defin ed a(±) := a(y'^^^). The set of 
Eqs.( ^.16 ) and (2.17) is exactly the same that was ob- 
tained in Paper I once we specialize the background bulk 
geometry to pure anti de Sitter. 

Using Eq. ( 2.17 ), it can be easily seen that the trace 
part of the metric junction condition is trivially satisfied. 

The remaining junction condition is the one for the 
scalar field, Eq. (2.£). After some computations, it re- 
duces in the "Newton gauge" to 

3 ^ ' a^'ipo 

(2.18) 



where we have defined 



y"(±) T 2((^o/^o) 



(2.19) 



Combining Eq. (2.18) with Eq. (2.7), we obtain the equa- 
tion for 6^ , 



(7=± 



□ W^JV ^1 + XI 2e^''^'^(2/ - y^"^)^ ■ (2.20) 



III. MASS SPECTRUM 



A. Zero modes 



Let us first consider the solution of the source free 
equations by setting T^^"* = 0. We will first consider 
the zero eigenmodes of which are the most impor- 
tant because they correspond to the 4-dimensional mass- 
less field responsible for the propagation of a long range 
force. 

If we set □(4) = in Eq. ( |2l^ ), the solution for h'^FJ'' 

is 

hZ^^^h^^}ixna'iy) + h(^}{x^)a\y) f J^y (3.1) 

where h^},} and h^j^} are 4-dimensional TT tensors inde- 
pendent of y satisfying □ '^4) = q , Furthermore, the 
jimction condition ( 2.14 ) gives 



"■/it/ 



-2a 



(± 



(3.2) 



Hence, there is no extra constraint on h[}J , but h[fj must 
be written as the second derivative of a scalar function. 
The former mode corresponds to 4-dimensional gravita- 
tional wave perturbations. Although there seem to exist 
5 independent degrees of freedom in this mode, three 
of them are pure gauge. On the other hand, the latter 
mode should be classified as scalar- type perturbations. If 
we forget about the stabilization mechanism, this mode 
corresponds to the mode called radion in Ref. 12^ . 



Notice that Eq. (3^) gives a relation between ^1"^^ and 



^5 



^5 _ 2 



(3.3) 



In the present case with stabilization, ^^_|_^ also appears in 
the junction condition for the scalar-type perturbations. 
Hence, the scalar field perturbation must also be chosen 
to be compatible with this condition. Here we should 
note that this relation holds only when we restrict our 
considerations to zero modes. In the general case, which 
will be considered later, there is no reason for such a 
relation to be satisfied. 

Let us now consider the solution of the zero eigenvalue 
scalar-type perturbation. One solution is 



3 



This mode satisfies the scalar field junction condition 



( pl8|) with Q if /(I) = -2a^±)^(5j^). However, this 



solution can be transformed to nothing by a gauge 

transformation with parameters = a~'^f'^^^ /2, ^'^ = 

-/.i'V'«"'(y')7'''^(y')'^y'- 

The other solution is 



(3.5) 



How ever, this mode is not compatible with condition 
(3.3). Hence, if we include the stabilization mechanism, 
no physical massless mode is present in the scalar-type 
perturbation spectrum, in contrast with the the case dis- 
cussed in Paper I. 



B. Massive modes 

Let us now consider the lowest massive mode. They are 
also important because they give the leading order cor- 
rection to the zero mode truncation approximation. Fur- 
thermore, if tachyonic modes are present, such a model 
must be rejected. 

To estimate the lowest mass eigenvalue in the general 
case, one needs numerical calculations. In order to obtain 
analytical approximations, here we assume that the effect 
of the bulk scalar field back reaction to the background 
geometry is small. Namely, we assume 



3i/2 



< 1. 



(3.6) 



For the metric, we can use the pure anti-de Sitter form 



a{y) = 

For simplicity, here we set 



y 



(+) - 



0, y'^-^ = d. 



(3.7) 



(3.8) 



First we consider the tensor-type perturbations. The 
expression for the mode functions Ui{y) for the operator 
L'-^-^) satisfying the junction condition [dy — 2H)ui(y) = 
on the positive tension brane is found in Refs. |l4|,|ll|. 
Denoting the eigenvalue corresponding to Ui{y) by m~, 
the mode function is given in terms of Bessel functions, 
u,{y) cx {Ji{m,t)Y2{m,iey/') - Y^{m,£) .Um^ey ' ')} . 
The boundary condition on the negative tension brane 
reduces to 

{ Ji {md) Yi (mM/'^ - Yi {m,£) Ji (mM^'^ } = 0, 



and determines the discrete eigenvalues m^. For small 
rrii, this condition becomes Ji{miie'^/^) « 0, and hence 
(3.4) the eigenvalues are given by « e~'^^^jii~^, where ji is 
the i-th zero-point of Ji . Hence the physical mass of the 
lowest massive KK mode on the positive tension brane is 
given by « 3.8e"''/^f "1 while that on the negative tension 
brane by « 3.8£~^. 

Next, let us consider the scalar-type perturbations. 
When we discuss the scalar- type perturbations, we 
should not assume an scale factor of the form (3.7) 
from the beginning, because the non-vanishing second 
or higher derivative of a{ y) can be important. Hence we 
will derive formula (3T9) for the lowest mass eigenvalue 
without assuming ( |3.7D . Only to obtain an estimate of 
this expression, we will use the scale factor (3.7). 

To discuss scalar- type perturbations, it is convenient 
to introduce a new variable defined by 



q ■= 



2k a' 



(3.10) 



where a prime ' denotes derivative with respect to the 
conformal coordinate z defined by dz = a{y)~^dy, and 
we have introduced A = ±a^/^(/5Q. The signature in the 
definition of A is chosen so that A is continuous on the 
brane s. The junction condition for q is obtained from 



2[d. + iA'/A)]q^±2Aif^^ {y^y^^^±). (3.11) 
The perturbation equation in terms of q becomes 



1\" 2k ,2 



A) 3 
Y^±2Ai\^^5{y~y^^^). (3.12) 



Integrating the above equation once at the vicinity of 
the branes, we correctly reproduce the junction condi- 



tion ( |3.11| ). 

To obtain the solution of Eq. ( [3.12 ) , we will construct 
the Green function for the differential operator appearing 
in it. It is given by 



(2^)4 ^ m^ + k' ' 



(3.13) 



where is the mass eigenvalue. The mode functions 
satisfy 



A 



2k 



'■mjqi, 



(3.14) 



with the boundary condition [dz + (^'/A)] qi — 0, at 
In the weak back reaction case, the term 5U :— 



(3.9) 2K(y9Q /3 in Eq. (3.14) is assumed to be small. If we 



4 



perturbatively expand the mode function and the mass 

the lowest eigenmode 
for the unperturbed system is trivially given by 



eigenvalue with respect to powers of SU like qi — 



(1) , A (0) , (1) , 

q: +■ ■ ■ and rrii = m.\ + 



(0) 

% 



A'- 



(3.15) 



with {itSq^)'^ — 0, where A/" is a normalization constant. 

It is straightforward to obtain the next order correc- 
tion. We find 



A 



A^dz' 



dz" 



2k 



(3.16) 



In terms of Sip, the above first order correction is given 

by 



5ip^^^ =Afipo{y) / dy' 



(m«)2 



(3.17) 



which becomes negative. Although we cannot give a 
proof here, this seems to be the case in general. Here 
we simply conjecture it. If tachyonic modes appear, they 
mean the breakdown of the model. Hence, we do not 
consider the case in which ipQ vanishes at some point, 
and hereafter we assume that ipo has a definite sign. 

Without explicitly specifying any model, we can make 
a crude estimate of the lowest mass eigenvalue. To eval- 
uate the nume rato r in the r.h.s of Eq. ( 3.19| ), we can 
use the form (3.7) for the scale factor a. To evalu- 
ate the denominator, we use the following fact. Be- 
sides some exceptional cases, a^(^Q will have a minimum 
(cases with no minima will be discussed later). There 
</3o + 2H(po = Vg — 2Hipo = 0, and accordingly the 
integrand has a rather sharp peak. We denote this 
point by j/c- At this point, the second derivative of 

Hence, 



'•'^ipQ is evaluated as (8i?^ 



2V^ - AH)a^ipl. 



we can approximate the integral like J{l/a*ipQ)dy 
[4iJVa4(F^)2]^^^Jexp[-(4i/2 + y^').=y.(2/ " yc?]dy. 
Under these approximations, the formula for the mass of 
the lowest eigenmode is obtained as 



The junction condition on the positive tension brane is 
already imposed due to the choice of the integration con- 
stant. Here, for simplicity, we have taken the e'^' — s- 
limit. 

This expression is seen to be regular even when (pQ van- 
ishes at some points. Denoting the point of vanishing ipo 
by yo, we can show that a^ipQ can be expanded around 
this point as 



60F 



(3.20) 



-I v=ya 



To proceed further, we consider the specific model for 
the bulk potential Vb{(p) discussed in Ref. |l8[, 



(3.21) 



4 • 2 

a ipg 



'v=yo 



yo) 



+a^{y-yo)^ + --- 



(3.18) 



For this model, the background solution for the bulk 
scalar field in the weak back reaction case is already given 
by 



where := i(4iJ^ + VB)ly=ao- From this, it is easy to 
see that the expression for Si£^ is regular at y = ■ 

The junction condition (3.11) on the negative tension 
brane gives the formula for the lowest mass eigenvalue. 



ml 



2n 


" r''' dy 




' ' dy ' 


T 


Jy{+) 


1 


JyW a^vl 



(3.19) 



This expression is positive definite if there is no point 
at which ipQ vanishes. However, when (p^ vanishes, the 
integral in the denominator changes its signature at 
that poin t. Hence, TOq can be negative. As seen from 



Eq. ( 3.18 ), the integrand does not have residue at this 
point. Hence the integral is independent of the way of 
modification of the integration path. It is expected that 
the dominant contribution to this integral comes from 
the region near the poin t y^, where Lp vanishes. By sub- 
stituting equation (3.18), we can approximately evaluate 
the denominator in Eq. (3.19) as 



/ 



dy 



{y-yoY + oP-{.y-yoY 



(3.22) 



where i^i = 2^-^+^/W^^T^P , = 2i-^ ~ ^J M~'^ + 
and 



Si 



'^(+) - '^(-)e 



(3.23) 



and are the values of (po on the positive and 
the negative tension branes, respectively. 

Let us briefly elaborate on the stabilization distance. 
For simplicity, we will assume that the ratio (^(_|_) 
is not extremely large, and also that all the input scales 
are similar, i.e. k ~ . The condition for weak back 
reaction is 



« 



1 



(3.24) 



Since Lp\ does not have its maximum in the bulk, it is suf- 
ficient to consider the conditions for weak back reaction 
on the boundaries. Then, we can evaluate the values of 
(yio on boundaries, 



5 



<^o(2/^"^) ~ VHP(-) - 4^^Vl + (Af2£V4)e''^''^(+). (3.25) In this 



The conditions for weak back reaction on both branes 
become 



2 / 4e'^2<i <^(^ 



4^-5/2 



1 



(3.26) 



The condition on the positive tension brane is satisfied 
by choosing {case A) <C with < or 

(case B) < with > £"3/2. 

In case A, the condition on the negative tension brane 
is automaticahy satisfied. In this case, the stabihzation 
distance is sensitive to the change of the values of the 
vacuum energy on the branes. Hence, without specifying 
the complete model, we cannot estimate the stabilization 
distance. If we consider the case in which M^f^ is not 
small, we find that it is natural to assume that the sig- 
nature of </'(+) is different from that of '/'(-)• If ^(+) and 
(/?(_) have the same signature, the ratio between them 
must be chosen to be extremely large to find a node-less 
solution of (po with sufficiently large brane separation d. 
In the case that ^P[+) and have different signature, 
the solution of ip^ becomes node-less for any choice of 
parameters. 

Case B is the case discussed in Ref. iQ. In this case, 
to satisfy the condition for weak back reaction on the 
negative tension brane. 



1 



-V2t 



In 



1 



(3.27) 



is required. Hence, 'P[+) and V'(-) must have the same 

signature. To realize a sufficiently large value of d/t, 

it is necessary that |</5(+)| is slightly larger than 

If we take the small IvP limit, the above expression for 

the stabilization distance coincides with that obtained in 

Ref. 

Let us now return to the computation of TOq. First 
we consider the c ase i n which a^(/3g has a minimum 
and the estimate ( 3.20| ) is valid. As we shall see be- 
low, case B is exceptional in this sense. From the 
condition that — 2Hipo ~ at y — yc, we have 
(,('^2-i^i)vc = V1B1/V2B2. Substituting this estimate of 



ijc into (3.2C), we obtain 



(-B1B2 



(3.28) 



To solve the hierarchy problem on the negative ten- 
sion brane, we need to set e'^''^ - {GeV) / W {GeY)= 
10^^. Hence, d/l k, 37. Therefore, it is rather 
natural to suppose that M'^id is larger than unity 
If |</3(_|_)| is not much larger than as as- 



sumed, we can approximate —B1B2 ~ 



case, the mass eigenvalue ttiq becomes 



cal squared mass on the negative tension brane becomes 
0(«;|(^(+)(^(_)|M2e-2^"''^(Vi+(*'^^^V4)-i)) ^hile that on 

the positive tension brane is the same as ttiq. As for 
|(/3(±)|, they must be smaller than £~^^'^ for the approxi- 
mation of weak back reaction to be valid. Also, there is 
a factor M^e-^^^'^'^Vi+f^^^^^-i) , which takes its max- 
imum value (0.2/£)2 at M « 0.33£~^. Hence, the mass 
scale on the negative tension brane tends to be smaller 
than the typical background energy scale £~^. 

Fi nally, we consider the special case in which a^(pQ in 



Eq. (3.19) has no minima. This happens when one of 
the terms in ipQ = I'lBie'^^^ + 1^2^26^'^*' can be totally 
neglected. This condition becomes 



(3.29) 



(3.30) 



In the former case we can set _Bi « 0, and in the lat- 
ter case we ca n set B2 ~ 0. Then, for the former case, 
formula ( 3.19| ) gives 



4k 



2 2 



-2d/t 



y/l + (Af2^2/4) 

_ -4^1 + (M2<?2/4)d/£ 



(3.31) 



It is easy to see that case B corresponds to this case. 
The model discussed in Sec. 4.1 of Ref. ||l^ with nega- 
tive small value of b also corresponds to this case (for the 
definition of b, see Ref. ||l^). The factor in the square 
brackets is almost unity. Substituting the approximation 
1^2 ~ —M^l/A, we recover the result obtained in [ p^ . 
Note that their and k are our 1/(4k) and £^^, re- 
spectively. 

For the latter case, we have 



4k 



,,2,„2 „~2d/l 



Vl-H(M2f2/4) 

4^1 + {K-Pe^/4)d/£ _ 



1 



(3.32) 



The model discussed in Sec. 4.1 of Ref. |ig|] with posi- 
tive small value of b corresponds to this case. To keep 
the mass sufficiently large, we have to choose ex- 
tremely large, which is not compatible with the weak 
back reaction condition. 



IV. RECOVERY OF EINSTEIN GRAVITY 

As we have shown in the preceding section, in two- 
brane models with stabilization mechanism, a physical 
massless mode is absent in the scalar-type perturbation 
spectrum, which is different from what happens in the 



6 



case without stabilization mechanism discussed in Paper 
I. This fact indicates that the resulting 4-dimensional ef- 
fective gravity can resemble Einstein gravity at linear or- 
der. Let us show how it can be recovered. 

To compute the induced metric on the branes, it is 
convenient to transform back to Gaussian coordinates. 



Thus, applying "minus" the gauge transformation (2.12) 
the induced metric on each brane is given by 



"■fj.u 



+ (4.1) 



Here we have decomposed h)i,u into two parts, the zero 
mode contribution and the KK contribution h^i^^\ 
The last term represents the residual 4-dimensional gauge 
transformation. As we have already noted, Eq. ( 2.16| ) 

and Eq. ( ^.17 ) determining h^ffj'\y^^'') and are es- 



sentially the same as the ones obtained in Paper I for 
the case without stabilization mechanism. There, it was 
found that the induced gravity approximated by the zero 
mode truncation becomes of the Brans-Dicke type . 

Here we only quote the result for (y*-^') up to 4- 
dimensional gauge transformation, 



□ (4) 



167rG( 



cr=± 



T 



where SttG'^-' := kNgJ^,^ and 



N 



a dy 



(4.3) 



Hence, the only possible source of an additional contribu- 
tion is that from h^f^^^ or 0^(i/'^)). Although we have 
just shown that there are no massless degrees of freedom 
in the scalar-type perturbations, we will find that the 
contribution from (/'^(y'-^') gives the correct long range 
force required to recover Einstein gravity. 

First we consider the weak back reaction case discussed 
in Sec. HI-B. From the normalization condition of (/q"' we 



obtain 



formula for the mass ( 3.19| ) 



2j,'dz{l/Ay 



. Comparing it with the 
we find a simple relation 



dy 

n2 



(e 



2d/e 



(4.4) 



Then, taking into account the contribution from the 
modes with the lowest mass eigenvalue alone, the Green 
function for q, Eq. ( 3.13D , is approximated by 



G„ 



d'^k 



3e 



^lA-\y)A-\y') 



(27r)4 2K£(e2rf/^ 
3A-\y)A~Hy') 



d'^k k'^e^^'^^"'" 



/c2 



(4.5) 



Using this Green function, and with the aid of (3.10) and 
(|3.12|), we find that 0^ is given by 



2t 



(-f 



(+) 



6 



[ml 



□ (4)]-i(a 



(+) 



t(+) + 



(4.6) 



Only the first term inside the square brackets gives a long 
range contribution to the induced metric. The propaga- 
tion of this force is essentially due to Cf±)- The source for 
the second term exactly becomes the trace of the energy 
momentum tensor of the ordinary matter field. Hence, 
the force due to this term becomes short-ranged with 
typical length scale ~ am^^. For the first term, using 
Eq. (2.17) for S,^j^y we finally obtain the contribution to 



the long range component of /i^^'' coming from (j)^ as 



"□(4) 



(4.7) 



(4.2) 

where we have used the fact that SttG*^^ 



/(I — e 2ci/^^ jj^ ^jjg ^ffQ^ back reaction case. 



Subs titut ing this into Eq. (4.1), and using Eq. (4.2) and 
Eq. ( ^.17 ), the zero mode truncation reproduces the for- 
mula for the linearized Einstein gravity, 



□ (4) 



CT = ± 



T 



(4.8) 



In the above derivation of Eq. (4.7), we have used sev- 
eral approximations. This derivation has the merit of 
having a result with a rather easy intuitive interpreta- 
tion. The gravitational field is propagated through the 
massless field ^ point far from the source T^^^ 

"(5) 

^(±) generates a cloud of metric perturbations through 
the interaction with the massive KK modes. However, 
the above approximate derivation is not completely sat- 
isfactory because some aspects of general relativity are al- 
ready tested with very high precision. Hence, we present 
an alternative and more complete treatment below. 

First we consider t o a pply the same trick used in 
the second line of Eq.(4^) to the complete Green func- 
tion containing all massive modes. Then the long 
range part of the Green function reduces to Gg ~ 
—6^{Ax^)^^qi{z)qi{z')/mf. This replacement is ex- 
actly valid when we focus on the long range force. Now 
one can notice that to use this Green function neglect- 
ing the terms corresponding to the short range force is 
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equivalent to solve the equation for (f>^ ( 2.20D by setting 
□ =0 from the beginning. 

For the case with = 0, we h ave already obtained 
t he general solutions to Eq. ( ^.20|) , i.e., Eqs.(3^) and 
(3.5). For convenience, here we quote the previous re- 
sults in a slightly different notation, namely 



with 



u± 



E 

<T = ± 



2H 



a''{y')dy'. 



(4.9) 



(4.10) 



Then the junction condition ( 2.18 ) with U^^"^ — Q deter 
mines /^^-' as /'^ 

back into the expression for (f)^ , we obtain 



T'^J^a^(±)£.l±)- Substituting these 



0^(y(±)) 



-2ffe7 



(±) 



2N 



(4-11 



cr = ± 



Adding this contribution to the contributions coming 
from the zero mode TT part and in Eq. (4J) for 
the induced metric on the branes, Einstein gravity at the 
linear order is recovered. 

It will be illustrative to give the next order correc- 
tion. The source term for the next order correction 
is -□(4),/,^ (1 + Y.a=± 2e*'"'^(2/ - y^"^)), where we have 
denoted the solution of the lowest order approximation 



( pll ) by (j)^ . By using Eq. ( pT7|) , □(4)0^ is evaluated 



cr=± 



u,{y)al)T^^\ 



(4.12) 



Then, by using the standard Green function method, 
we obtain 



<(!/) 



dy' - K 



y(-) ^0 \ 



, (4.13) 



where and K- are integration constants determined 
from the junction condition as 



K4 



i(±) 



1 ± 



H 



(±) 



(4.14) 



For simplicity, we take again the e^^'' limit. Then, 
we have K± = 0, and we obtain 



(4.15) 



where we have defined 



Ut{y)uj{y) 



dy. 



(4.16) 



In the weak back reaction case, we can approximately 
evaluate this integral by using the following facts. First, 
we note that the combination a'^u± can be shown to be a 
slowly changing function for y ^ yc 'm the weak back reac- 
tion case. As before, j/c is the peak location of (a'*(/3o)~^. 
To show this, we use the fact that ipo is approximated 
by viBie^^y for y ^ yc- This c onsta ncy of a^u± indi- 
cates that the integrand of Eq. ( 4.16 ) is approximately 
proportional to (a^(^§)^^. Then by us ing the background 
geometry defined by Eqs. (3/7) and (|3.8|), /+,+ is evalu- 
ated like 



\y=d 



dy 




where we have used (3.19). The other components are 
also evaluated in a similar way to obtain the relations 
/+,+ « e-2'^/^/+,_ « e-'^'^/^I--. The correction ob- 
tained by substituting these estimates for /^j into ( 4.15| ) 
is consistent with the contribution from the second term 
in the square brackets in Eq. (4.6). If we set D^^^ — 
there, we recover the same result. 

It is easy to check that the treatment of taking D^**^ to 
be small is consistent if the distribution of the energy mo- 
mentum tensor is sufficiently smooth. In the treatment 
presented here, we have taken into account the contri- 
bution from all the massive modes simultaneously. The 
approximation is completely valid as long as we consider 
smooth matter distribution as compared with the mass 
scale of the lowest massive mode. However, the lowest 
mass on the positive tension brane becomes very small if 
we consider the case in which we are living on the nega- 
tive tension brane. In this sense, the treatment presented 
here is rather restrictive when we discuss perturbations 
caused by the matter fields on the positive tension brane. 



V. TT PART REVISITED 

Applying the same technique that we have used in the 
preceding section for the scalar-type perturbations, we 
can also deal with the KK mode contribution for the TT 
part /i/jt/ . By using the Green function method, we can 
easily calculate the zero mode contribution like 



e^ = -2W(y)(nW)-i5]a^,)I](X), 



(5.1) 



cr=± 



tutmg ft^y 



After using Eq. ( 2.17 ), we can recover Eq. (12). Substi- 
+ h'^Jfy^^ into ( ^.16| ), we obtain the 



equation for the KK contribution, 



cr = ± 



^□(4) ^ liTT) 



hiKK) 
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= 2nY^ E^^) [al^N - S{y - . (5.2) 

As before, neglecting the □'^''^-term for the massive KK 
contribution, we can solve this equation like 

I'' dyV(y")-C(.)V (5.3) 




where C(+) and C(_) are constants. We should recall that 
we have already subtracted the zero mode contribution. 
Hence h^J^^^ must be orthogonal to the zero mode. This 
requires 

„(-) , „„(-) 



V 

(+) a 



dyhf,''\y)=0. 



(5.4) 



The constants C(_|_) and C'(_) in the solution ( pT^ ) are 

determined by imposing this cond ition . 

For simplicity we again adopt (3.7) with ( p.Sl) as the 
background geometry. Then, from condition ( |5.4D , C(j-) 
are explicitly calculated, and the resulting KK contribu- 
tion becomes 



4(1 



2d/fj L. 



a=± 



1 



- 1 



(5.5) 



Hence, on the respective branes, the contributions from 
the KK modes become 



dKK) 



((3-4d/^)EW+E(;) 



(5.6) 



where we have assumed djt^ 1. Here we should recall 
the remaining degrees of freedom for the 4-dimensional 
gauge transformation. Using these degrees of freedom, 
Ey^jy ca n be replaced with 

T/i with the aid 



of Eq. ( 2.11 ). Hence, the KK modes, of course, do not 
give any long range force contribution. If one takes the 
d/^ — > cxD limit, seems to diverge. But this is just 

due to the breakdown of the approximation. In this limit, 
the mass difference of the KK modes becomes zero. 



VI. SUMMARY 

In this paper we have developed a systematic proce- 
dure to evaluate the perturbations in the 5-dimensional 
brane world model proposed by Randall and Sundrum 



g supplemented with the moduli stabilization mecha- 
nism by Goldberger and Wise |l^ ]. 

We have first investigated in detail the mass spectrum 
of this model. In the case without stabilization mecha- 
nism, there was a scalar-type massless mode, which was 
called radion in Ref. |2^]. We have clarified how this 
massless mode disappears once we switch on the stabi- 
lization mechanism. 

We have also estimated the mass eigenvalue and the 
mode function corresponding to the lowest mass eigen- 
modes for both tensor-type and scalar-type perturba- 
tions, assuming that the back reaction to the background 
geometry due to the bulk scalar field introduced for the 
moduli stabilization is not large. The physical mass of 
the lowest tensor- type mode becomes ^ l^^ on the nega- 
tive tension brane and ~ e~'^^^l~^ on the positive tension 
brane. Here t is the curvature scale of the background 
geometry, and d is the proper distance between the two 
branes. In the original model, is supposed to be TeV 
scale. 

For the physical mass of the lowest scalar-t ype m ode, 
we have obtained rather general formulas, ( 3.19|) and 
( |3.20D . To proceed further in estimating the mass, we 
have specified a model for the bulk potential of the 
scalar field. We have considered a simple quadratic 
potential whose mass is given by M, and we have as- 
sumed that the 5-dimcnsional gravitational constant is 
also ~ . We pointed out that in this model there are 
two regimes in which the weak back reaction condition 
holds. The first case is the one in which the vacuum ex- 
pectation values of the scalar field on the branes, <^(±), 
are sufficiently small compared with the background en- 
ergy scale, i.e., '■P(±) ^ £~^^^. In this case, we found 
that the mass of the lowest scalar-type mode becomes 
VI'^(-)<^(+)|^^Me-(Vi+(A^^^V4)~i)<i/£. (For a more 



precise formula, see Eq. ( 3.28 ).) If we take into account 
the fact that we need to set d/^ « 37 to solve the hier- 
archy problem on the negative tension brane, this factor 
is at most « 0.2£-iy'|(^(_)(/3(+) for Mi w 0.33. The 
second case is the one in which M is small compared with 
i"^, which is the situation discussed in Ref. In this 
case, we reproduced the same result found there. Namely, 
the mass is approximately given by {\(p(^_)\i"^^'^)M'^£. As 
pointed out in Ref. ||l8|, the mass scale on the negative 
tension brane tends to be smaller than the typical back- 
ground energy scale £~^. We confirmed that this is a 
general feature within the context of weak back reaction. 
However, looking at these formulas, it seems that we can 
raise the mass by taking a slightly larger vacuum ex- 
pectation values of the bulk scalar field on both branes, 
although a large value of |v'(±)| is inconsistent with our 
approximation. Hence, if we remove the technical limi- 
tation of weak back reaction, it is not clear whether the 
physical mass of the lowest scalar-type mode is always 
smaller than that of the lowest tensor-type mode. We 
also mention that the mass scale on the positive tension 
brane is smaller by a factor of e^'^l^ . 
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Next, we have developed a method to evaluate the ex- 
plicit form of the perturbations caused by the matter 
fields confined on the branes. This is a generalization 
of the results presented in the previous paper jl^. We 
found that, as is expected, Einstein gravity is exactly re- 
covered for the long range force at the order of linear per- 
turbations when we impose the stabilization mechanism. 
The formulas for the leading correction from the scalar- 
type perturbations ( 1.15| ) and that from the KK modes 
(5.6) are also derived without specifying the model for 
the potential of the bulk scalar field. From these formu- 
las, we can read the coupling of the metric perturbations 
induced on the branes to the matter fields on both branes. 



As we have confirmed in Sec.V, there is no pathological 
behavior in all perturbation modes at the level of linear 
perturbations. This is a rather expected result because 
almost all important information at the level of linear 
perturbations is contained in the mass spectrum except 
for the strength of coupling to the matter fields. How- 
ever, it is still unclear what happens once we take into 
account the non-linearity of gravity. To investigate this 
issue, the formulas obtained in this paper will be useful. 
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